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ABSTRACT 

We discuss the structure of clusters in a class of flat cosmological models with the 
fraction of mass ftcDM — 0.8 in cold dark matter, and the rest in hot dark matter in the 
form of massive neutrinos. These Cold + Hot Dark Matter (CHDM) models with one, 
two or three massive neutrinos, with total mass ~ 4.6 eV, are considered. Neutrinos 
of that mass range are too hot and cannot constitute the halos of galaxies and groups, 
but only of clusters of galaxies. The limit on the density of neutrinos in the central 
parts of galaxy clusters is estimated from the phase space density constraints. The 
ratio of the density of neutrinos to that of cold dark matter through the cluster is found 
analytically. It appears that the density of the neutrinos is suppressed within the Abell 
radius. However, neutrinos contribute ~ 20% of the mass density to the cluster halo. 

Our numerical simulations fairly match analytical results. The simulations indicate 
that the cluster halo dark matter density profile has the power-law slope a ~ 2.5, which 
is close to that in the model with the cosmological constant. We also found that in 
the CHDM models the velocity dispersion is almost constant across the cluster. This is 
quite different from the model with the cosmological constant or from the open model 
where the velocity dispersion in the cluster outskirts declines. 

We discuss observational tests that can probe cold, neutrino, and baryonic compo- 
nents unequally distributed in clusters: X-ray emission and weak gravitational lensing. 
We input the found spherically symmetric fit to the CHDM mass density profile and 
the X-ray surface brightness for the cluster A2256 into simple model of the hydrostatic 
equilibrium of the hot gas. The derived theoretical temperature around the center de- 
partures from both the observational data and actual prediction of the cosmological 
model, which give almost constant temperature. We also find that the problem of high 
baryonic fraction in clusters is not resolved in the CHDM models. 

Subject headings: theory - dark matter- clusters of galaxies 
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1 Introduction 



Most cosmological models are based at present on variations of the Cold Dark Matter scenario, 
with initial conditions compatible with inflationary theory. Two currently popular models of that 
type, which satisfy a range of observational data, are the Cold + Hot Dark Matter (CHDM) 
model and Cold Dark Matter model with a cosmological constant (ACDM). These models have 
an extra parameter, f — ilcDM^ which can be adjusted to make the models compatible with the 
power spectrum of density fluctuations P{k) estimated from galaxy catalogs, large-scale velocity 
streaming, and CMB AT/T fluctuations. Another set of constraints comes from the abundance 
of protogalaxies at high redshifts and clusters at the present epoch. The CHDM model assumes a 
fraction (20%-30%) of the critical density to be in the form of massive neutrinos. In the ACDM 
model most of the mass of the Universe (60%-80%) is due to the cosmological constant A. In these 
respects the models are relatively well studied, so that one can put limits on the region of the 
parameter space in the (I — ilcDM^h) plane, h being the Hubble parameter in units of 100 Mpc 
km~^ s~^. An earlier version of the post COBE CHDM model adopted Oj, = 0.3 from one massive 
(supposedly r-) neutrino with mass = 7 eV ( Davis et al. 1992; Klypin et al. 1993). The allowed 
region of the parameter space (Oj,, h) was studied by Pogosyan & Starobinsky (1993) and Liddle & 
Lyth (1993); protogalaxy abundances were studied by Mo & Miralda-Escude (1994); Kauffman & 
Chariot (1994); Ma & Bertschinger (1994); Klypin et al. (1995), and clusters abundance and their 
statistical properties by Klypin & Rhee (1994); Jing & Fang (1994); Croft & Efstathiou (1994); 
Cen k Ostriker (1994); Bryan et al. (1994); Walter k Klypin (1995). At present, the model with 
values (Oj,,/i) = (0.2, 0.5) has been adopted. A disadvantage of this model is that it is in the 
conflict with the high value of the Hubble constant (e.g., Ereedman et al. 1994). 

Some interpretations of current neutrino experiments do not exclude the models of two massive 
neutrinos with close masses ~ 2.3 eV each (Primack et al. 1995). Erom the general point of view, it 
is instructive to consider even more exotic combinations with three massive neutrinos. The density 
parameter of the Universe associated with the massive neutrinos is Oj, O.OII/i"^ (?7ij/,/IeV). 
Because the neutrino Jeans length depends on the neutrino mass (Doroshkevich et al. 1980; Bond 
k Szalay 1983), the shape of P{k) in the CHDM models depends not only on the total neutrino 
density Oj,, but also on the masses of neutrinos of different types. Thus for fixed total density Oj, 
the spectrum P{k) will depend on how the mass is distributed among three types of neutrinos. 
Here we will consider the general range of neutrino masses, but under the constraint that its total 
value is Y^ii^iy, = 4.6 eV, which corresponds to Oj, = 0.2 (for h = 0.5). In the CHDM models with 
slightly tilted primordial spectra (Pogosyan k Starobinsky 1995a) or with more than one massive 
neutrino (Pogosyan k Starobinsky 1995b), the range of [l — ilcDM^ h) parameters is slightly broader 
than for the standard CHDM model with one neutrino. However, the allowed value for the Hubble 
constant is h < 0.6, which still might be too small. Note that for the ACDM model the allowed 
parameters are given by ilh 0.25 and > 0.3 (Kofman, Gnedin k Bahcall 1993; Ostriker k 
Steinhardt 1995). 

To study the viability of the CHDM model against the accumulating observations of clusters, 
such as X-ray emission, strong lensing producing arcs and weak lensing producing background 
galaxies distortion, and galaxy velocity dispersion, one needs to know the mass and velocity dis- 
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tributions in the dark matter components in the clusters. In this paper we study the density and 
velocity profiles of clusters in the CHDM cosmological models. We will make calculations primarily 
for the CHDM models with one 4.6 eV or two 2.3 eV neutrinos. However, we would like to note 
that the physical effects we are discussing here are more general and are applicable to a broad range 
of the CHDM-type models. We also compare the results with the corresponding predictions in the 
ACDM model, and make preliminary comparisons with some observations of clusters. 

It is well known that the properties of clusters such as the mass density profile p[r) and the 
(galaxy) velocity dispersion profile cry[r) can probe the background cosmological model, as well 
as the index n of the power spectrum P{k) (see, for example, Gunn 1972; Hoffman & Shaham 
1985; Quinn, Salmon, & Zurek 1986; Crone, Evrard, & Richstone 1994 and references therein). For 
hierarchical models like CDM, the formation of a cluster can be treated, at least qualitatively, as a 
spherical accretion onto a peak of the initial density field. For the peak one can take, say, a massive 
galaxy with mass IO^"^M0 (Gunn 1972). The accretion model gives us insight into how the density 
and velocity profiles in the cluster outskirts depend on the parameters of the cosmological model 
like the CDM model or its interesting variations. 

However, in the CHDM model the process of clusters formation is more complicated due to 
the nature of the neutrino component. Indeed, hot neutrinos accrete onto the protocluster only 
after their thermal velocities drop below the velocity dispersion of the protocluster. Thus, the 
neutrino halo of the cluster starts to form after the cold dark matter halo is well advanced into the 
nonlinear stage. Neutrinos also cannot be compressed to an arbitrarily high density in the cluster 
center. Because neutrinos are coUisionless particles, their phase-space density is conserved. This 
puts a limit on their mass density in gravitationally bound objects. The sizes of the objects that 
should have the neutrino halo are related to the neutrino masses. For the pure hot dark matter 
model the phase space density constraint leads to the Tremaine & Gunn (1979) criterion, which 
links the mass of the neutrino to the parameters of the isothermal halo: the 24 eV electron neutrinos 
compose the halo of objects of the size of normal galaxies. One expects that in the CHDM model, 
a few eV neutrinos do not contribute to masses of galaxy halos, but they can make a contribution 
to masses of clusters. In this paper we extend the Tremaine-Gunn criterion for the CHDM models 
and will apply it for the clusters of galaxies. 

In this paper we also address two specific questions of the cluster structure for the CHDM-class 
models. How small is the fraction of the neutrino mass density, / = p^/pcDM^ in the core of a 
cluster? This question is rather topical, because it is related to the problem of the baryon-to-dark- 
matter ratio in clusters. We calculate the compression factor / as function of the neutrino masses. 
Next, what is the profile of the neutrino mass density in the halo of clusters relative to that of 
CDM? This question is related to other important problems: the potential discrepancy between 
the total surface mass density which can be reconstructed from weak gravitational lensing, and the 
X-ray emission from clusters. 

In Sec. 2 we study the differences in the evolution in the linear regime between the cold dark 
matter and the neutrinos. In Sec. 3 we derive the phase space density constraint on the neutrino 
density for virialized clusters. In Sec. 4 we derive the cold particles and neutrino density profiles in 
clusters, first analytically and then by N-body simulations. In Sec. 5 we discuss some observational 



2 



constraints on the baryonic, cold dark matter and neutrino components in clusters. First we derive 
the temperature profile of the hot gas in the hydrostatic equilibrium with the CHDM mass density 
and compare it with the data. Then we obtain the total surface density in the CHDM model which 
is probed by the weak gravitational lensing. Finally we discuss the galaxy velocity dispersion profile 
in clusters as a discriminative test for cosmological models. 

2 Linear Fluctuations in CDM and Neutrinos Compo- 
nents 

In a model with only a cold dark matter component, the formation of a protocluster can be roughly 
viewed as an infall of spherically symmetric shells onto a local maximum of the primordial field of 
gaussian density perturbations 6[r). In order to study some statistical aspects of cluster formation, 
such as cluster abundance with redshift, the primordial density field is smoothed with a filter, with 
a scale corresponding to the cluster mass, ~ lO^^M©. Then the cluster abundance can be estimated 
via the statistics of the overdense regions of the filtered (5-field above some threshold (5c, as in the 
Press-Schechter model. In this paper we address instead the shape of density profile of clusters. 
For this purpose we consider the initial density peaks relevant for the protoclusters, and smooth 
the initial density field on that scale, which corresponds to a massive galaxy with mass 1O^"^M0. 
The structure of a cluster halo forming by the accretion onto the initial overdensity is then defined 
by the details of the density profile around the maximum. The later can be well described by the 
two-point linear correlation function ^(r) =< 6[r)6[0) > (e.g., Hoffman 1988). 

Specific features of cluster formation in the CHDM model with cold and hot components can 
be calculated in the linear regime. Let us consider the density profiles in both components around 
the density peak. Linear evolution of the adiabatic perturbations in the CHDM model is derived 
from the kinetic equation for neutrinos and the linearized Finstein equations (see Appendix A). 
Solving these equations numerically, one can calculate the transfer function of the total density 
fiuctuations, which is C{k) = k~^P[k) for the scale-free initial power spectrum. Figure 1 shows 
transfer functions C{k) at 2; = for three CHDM models: only one neutrino is massive, = 4.6 
eV (upper curve); two neutrinos have equal masses 2.3 eV (middle curve), and all three neutrinos 
have equal masses 1.53 eV (lower curve). The upper and the lower curves are the limiting cases. 
Curves corresponding to all other distributions of 4.6 eV among three types of neutrinos will be 
sandwiched between these limits. This rule is valid not only for the transfer function but for other 
characteristics we will consider throughout the paper (see also Pogosyan & Starobinsky 1995b). 
The splitting of C{k) for different masses on clusters scales is quite noticeable. 

The overdensity profiles in cold (A = 1) and hot [A = 2) components around the peak can 
be expressed through the corresponding correlation functions, (5^(r) = z/(^^(r)/(T^, where p is the 
peak height in units of the rms density dispersion a a- The density profile is simply /Oa(?") = 
0^/j(l -|- (5^(r)). From the linear equations of Appendix A we numerically calculate the two-point 
correlation functions for both components, as functions of the filtering mass M and redshift z. For 
definiteness, we choose the model with two 2.3 eV neutrinos. In Figure 2 we plot average CDM 
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(solid curve) and 2.3 eV neutrino (dashed curved) density profiles /Oa(?") around a density peak; 
in upper panels for filtering on the scale fO^"^M0 (with the peak heigh v = f.68), in lower panels 
after filtering on the scale fO^^M©; in the left panels at 2; = 4, in the right panels at at 2; = 0. The 
difference in the evolution of the overdensities at a given scale depends on how that scale relates 
to the neutrino Jeans length Aj: 



This corresponds to the Jeans mass Mj ~ fO^"^M0 at 2; = 4 (for 2.3 eV neutrinos). For the upper 
left panel in Figure 2 the scale of the fiuctuations is less than Aj. As the result, the growth of 
the CDM fiuctuations at that scale is significantly advanced relative to the neutrino fiuctuations. 
Compare this with the evolution shown in the lower panels at scales greater than Aj, where the 
growth in both components are comparable. Therefore the formation of clusters begins with the 
formation of the CDM protoclusters. 

At early stages hot neutrinos do not accrete on the CDM protoclusters because the typical 
neutrino velocity Vi, is too large: 



Note that Aj ~ / dt-v^. Neutrinos accrete on the protocluster only after drops below the velocity 
dispersion of the protocluster. At present, neutrinos have accreted into the halos of clusters. This 
problem is reminiscent to the neutrino accretion on the nonperturbative density inhomogeneities 
(seeds) formed by cosmic string loops (Brandenberger, Turok & Kaiser 1986). If one models a seed 
as a delta- function overdensity, then neutrinos form an extended halo on the scale ~ Aj around it. 

One might think that in the CHDM model, there would be an extended neutrino halo, where 
neutrinos dominate over cold particles, around a cluster where cold particles dominate over neutri- 
nos. This is not true. In the CHDM model neutrinos make a contribution to the cluster halo, but 
their density nowhere dominates over the CDM density. As we will see from N-body simulations in 
Sec. 4, the ratio of the neutrinos density to the CDM density never exceeds the background ratio 
1/4. This result can be understood if we consider the spherical infall model implemented for the 
neutrinos whose thermal velocity is already much lower than the clusters velocity dispersion. 



For clusters in equilibrium, a useful constraint on the density of the neutrinos inside massive halos 
can be derived from the distribution of neutrinos in phase space. In this and the next Sections we 
extend the Tremaine & Gunn (1979) criterion to the CHDM model. The constraint comes from 
the condition that for the coUisionless particles the phase space density n(p, x) is constant along 
the particles trajectories. The initial phase space density of neutrinos of one family Pi is given by: 







3 Neutrino Phase Space Density Constraint 




(3) 
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where is the temperature of neutrinos and p is their momentum. The initial phase space density 
distribution riinit has the maximum value ninit,ina.x = 0.5 at p = 0. Thus by equating the final phase 
space density with the maximum of the initial phase space density we get upper limit on the density 
in the final distribution. 



Let x) be the final phase space density of neutrinos of one family Pi gravitationally bound 

in the massive spherically symmetric halo. The space density of neutrino Pi at radius r is the 
integral of x) over the momentum p = m^^v 

p.Ar) = ^,fd'pnMr). (4) 

Let the one-dimensional velocity dispersion in the halo be cr, then in (4) we can factor out all 
dimensional parameters 

4 3 

where y = v/cr. Now we can estimate the maximum value of the neutrino Pi mass density at the 
center p^^[0). The final phase space density is assumed to depend only on the energy 

nf{p,x) = no ■ f{E), (6) 

where Uq is a normalization constant, and E is the particle energy: E = + cf){r) j a"^ . The 

gravitational potential cf){r) is normalized at the minimum to ^(0) = 0, and the distribution func- 
tion f{E) is normalized to /(O) = 1. The integral in (5) for common distributions f{E) (including 
those we will use below) can be evaluated as / dyy^ f[y^ /2) ~ 0(1); for example, for the Boltzmann 
distribution / dyy^ f[y^ /2) = 1.25. Below we will distinguish the fine- and coarse- grained distri- 
butions, see e.g. Binney & Tremaine (1987). The coarse-grained distribution < n[p^x) > is the 
fine-grained distribution averaged over some phase volume; when the pair relaxation is negligible 
that is always lower than the actual, fine-grained space density. 



The phase space density is invariant along the neutrino trajectories. Therefore the final fine- 
grained distribution (6) cannot exceed the maximum value 1/2 of the initial distribution (3), Uq < 
0.5. Thus, an upper limit on the neutrino Pi density in the halo center is 
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«,.(0)<^^. (7) 



We can give not only the upper limit but also a reasonable estimation of the central neutrino 
mass density. An average (coarse-grained) initial phase space density < n(p, x)init > can be found 
either from the mean hot neutrino energy, or from the median hot neutrino energy, or from the 
condition that half of all neutrinos exceeds the value < n[p^x)init >■ All these methods give us 
about the same value < n[p^x)init >~ 0.05. By equating the final average phase space density to 
this value, we obtain the expected neutrino density at the center 

4 S 
TTl (7 

.,,(0)^0.1.5^. (8) 
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Using the equation (8), we can evaluate the neutrino mass density at the center of clusters 
for different neutrino masses. In the model with a single massive neutrino with mass 4.6 eV, we 
nave p,y (0) ~ 2 • 10-26 (m^/4.6eV)^((7/fOOOA;m/s)''5/cm3; in the model with two massive neu- 
trino with mass 2.3 eV each, the total central neutrino mass density is /Oj/(0) = 2/Jj,^(0) ~ 2.5 • 
fO-2^ (m^/2.3eV) {a / l{){){)km / sf g / cm^ ■ in case of three massive neutrinos with mass 1.53 eV each, 
total central neutrino density is /Oj/(0) = 3/Jj,, (0) ~ 7.4 • IQ-^* (mj,/I.53eV)^ [a / 1000km / s)^ g / cm^ . 
The neutrino mass density at the center strongly depends on the neutrino mass. However, it is only 
by factor 10^ — 10"^ greater than the mean density p = 1.88 • 10~^^h^ g/cm"^. The CDM component, 
filtered with the scale of a big galaxy, has at the center pcDM — 3 • lO'^p. Thus, at the center, the 
CDM component significantly dominates over the HDM one. 



4 CDM and Neutrinos Halo Density Profiles 

We study the neutrino density profile Pu{r) together with the CDM density profile /Ocdm(?") in 
the clusters. An especially interesting characteristic is their ratio /(r) = Pi,{r) j pcdm{''^)^ which 
describes the relative compression of the cold particles in the clusters. First, we use an analytic 
model based on certain assumptions about the neutrinos and cold particles phase space distribu- 
tions. Then we compare our predictions with results of N-body simulations for the CHDM model. 



4.1 Analytical Model 

Finding the mass density profiles of the two components requires a generalization of the one (neu- 
trino) component problem considered by Tremaine & Gunn 1979. To derive the neutrino density 
profile Pi,{r) together with the CDM density profile Pcdm{''^) in the clusters, in the first turn, one 
should specify the phase space density distributions in both components. 

A convenient model for the CDM particles is the spherically symmetric isothermal distribution 

ncDMiv.r) = N ■ eyi^i--^], (9) 

V <old/ 

where E = v"^ l2 -\- cf>{r), Ccoid is the velocity dispersion of cold particles and is a normalization 
constant. For this phase space density distribution, the mass density distribution of the CDM 
particles is 

( 

PcDM[r) = Pcdm{0) ■ exp ( — ) , (10) 



Coldy 



where Pcdm{^) is the central CDM density. 



We take the following distribution function f{E) for neutrino (see eq. (6)): 

I + exp (L/a^ - p) 
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where cr^ is the neutrino velocity dispersion, and fi plays the role of the chemical potential, which 
we treat as a free parameter and will be fixed below. Note that the cold particles velocity dispersion 
Ccoid, generally speaking, is different from cTj,, but here we assume they are equal. Using (f f) in the 
formula (5), we obtain the neutrinos mass density p^^[r)^ see Appendix B for technical details. At 
this point the parameter fi can be fixed by considering the ratio Pi,{r) j pcdm{''^) at large distances 
r where it should match the background ratio Oj,^ : VIcdm- Then the neutrino density profile Pi,^{r) 
is the function of the gravitational potential cf){r) and a parameter /Ocdm(O) only, see eq. (27) of 
Appendix B. 

Now we are ready to write down the Poisson's equation for the gravitational potential cf){r) for 
the two component spherically symmetric configuration 

^ + -^ = 4vrG(^aDM(r) + /..(r)), (f2) 

where the CDM and neutrino densities as functions of cf){r) are given by eqs. (fO) and (27). In 
Appendix B we reduce the Poisson equation to the form of a generalized Emden equation for 
a selfgravitating isothermal distribution. We solve numerically the Emden equation for the two 
component halo for given parameter /Ocdm(O). We have to choose the value of /Ocdm(O). Generally 
speaking, how large can this value be is still an open question, since for the CDM models density 
has I/r cusp up to the scale of resolution (Dubinsky & Carlber I99I). Our numerical simulations 
have a resolution of a large elliptical galaxy. Therefore as a central CDM density in clusters we will 
mean the "coarse-grained" density. We choose /Ocdm(O) equal to 3 • 10'*/?, the average density of a 
large elliptical galaxy. 

In Figure 3 we plot the final results in terms of the ratio of the two density profiles /(r) = 
Pi,{r) I pcdm{''^)^ for three cases: one 4.6 eV, two 2.3 eV and three 1.54 eV neutrinos. We see that 
in the case of two 4.6 eV neutrinos the neutrino mass density in the core is suppressed by a factor 
of about 2 relative to the CDM density. In the case of one 2.3 eV neutrinos the suppression factor 
is about 5. In the case of three 1.54 eV neutrinos the suppression factor is about 10. At distance 
about I Mpc, neutrinos contribute to the cluster outskirts in the proportion Oj, : VIcdm- In the 
limit of large r, the isothermal analytical model, obviously, gives for both densities the asymptote 
oc r~^, which is, as we will see, is slightly fiatter than the actual dependence oc r~^'^, which we 
find in numerical simulations. However, the analytical model predicts the ratio of the two density 
profiles fairly accurately comparing with that from numerical simulations (upper panel of Fig. 4). 

4.2 Numerical Simulations 

To obtain detailed density structure of clusters in the CHDM model, we made an N-body simulation. 
In addition to the density profiles, we use the simulation to get the profiles of the velocities in 
the clusters. The simulation has been made in a box of 50/i~* Mpc size with periodic boundary 
conditions. It has 256"^ cold particles and twice that many hot particles moving in a 800"^ mesh. The 
(formal) resolution is 62.5/i~* kpc. Initial conditions correspond to a COBF normalized (Qrms-ps = 
I8/iK) CHDM power spectrum with total = 1, ^cdm = 0.8, Oj, = 0.2, and h = 0.5. We chose 
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the CHDM model with two equal mass 2.3 eV neutrinos. We also will compare the results with 
those for the "old" 7 eV model. For this case we use another simulation with the same number of 
particles and box size, but with slightly worse resolution (512"^ mesh). 

Two largest clusters - (a) and (b) - in the simulation for the model with two neutrinos were 
chosen as examples. They are smaller, but comparable to the Coma cluster. The cluster (a) 
has mass 6.5/i~^lO^'*M0 (left panel in Figure 4) and cluster (b) has mass 5.5/i~^lO^'*M0 (right 
panel in Figure 4) within the Abell radius 1.5/i~^Mpc. In the lower panel of Figure 4 we plot the 
corresponding CDM density profiles. Now we discuss the fit by simple analytic form for the CDM 
density profile. There are uncertainties in the fit due to the numerical resolution of our simulations: 
the very central part of the clusters is not resolved. 

Navarro, Frenk, & White (1995) found that for the CDM model the dark-matter density profile 
in central parts of clusters is well approximated by p[r) = 60/jr^/r(r + 0.2r^,)^, where Vy is the 
virialized region the cluster (radius at which the mean overdensity is equal to 200). We find that 
this formula also gives a very good fit to the CDM density profile in our CHDM models inside 
central lh~^ Mpc region. But at larger radii the density profile in CHDM models is shallower: 
it declines as oc r~", where a is usually less than 3. To fit the density profile both within and 
outside cluster, we will use a simple conjecture which simultaneously gives us 1/r asymptote at the 
center, and 1/r" asymptote on the cluster outskirt 



This density profile gives a better fit than the King-like profile, although the later is also quite 
reasonable. The profile (13) provides a central surface mass profile large enough to produce arcs 
from the strong gravitational lensing. Fquation (13) has been tested for several clusters in the 
simulation. We find the slope in the halo for the model with two 2.3 eV neutrinos to be a 2.5. 
The same slope 2.5 in the halo was earlier found for clusters in the CDM-|-A model (Crone, Fvrard 
& Richstone 1994). Thus, clusters CDM density profiles slopes in the 2.3 eV CHDM and CDM -^A 
models are similar. The value of the scale parameter R in (13) we found is i? ~ 200 — 300 kpc (for 
h = 0.5). The neutrino density profile is fitted by the formula /0„„(r) = f{r)pcDM{''^)^ where /(r) 
is the compression factor, plotted in Fig. 3. 

In Figure 5 we plot the CDM density profiles in the "old" CHDM model with nii, = 7 eV. 
In this model the density profile slope is steeper, a sa 2.15. It is interesting that the slope a is 
rather sensitive to the details of the CHDM models, such as Oj, and the spectrum of the neutrino 
masses. Partly the difference in a is caused by different power spectrum index n at the clusters 
scale. For Oj, = 0.2, 2.3 eV neutrinos, n sa —1.36; for Oj, = 0.3, 7 eV neutrinos, n sa —1.56. The 
upper panels of Figure 4, 5 show the ratio /(r) of hot-to-cold densities as a function of distance 
from the cluster center, normalized in such a way that unity corresponds to the cosmological ratio 
Oj, : Vic DM- We see that indeed the neutrino density at the center is suppressed by a factor 5 
relatively to the CDM density in the case of 2.3 eV neutrinos (Figure 4); and by factor 2 in the 
case of 7 eV neutrinos, in agreement with the analytic prediction, cf. Figure 3. At about ~ lh~^ 
Mpc the ratio /(r) saturates and reaches its background value Oj, : VIcdm- The fiuctuations of the 
ratio /(r) at large r are not statistical. Both clusters (a) and (b) have neighbors - smaller clusters 




(13) 
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at about Ih ^ Mpc from the cluster (a) and 2h ^ Mpc from cluster (b). As the result, the minima 
of /(r) correspond to the positions of other clusters and groups. 

What appears to be most interesting is the profile of the velocity dispersions Urms(?") in the 
clusters. Figure 6 shows Urms(?") for clusters (a) and (b) for the 2.3 eV model. The solid curve 
shows Vrmsi'i') of cold particlcs; the dashed curves are for hot particles. Outside of the Abell radius 
both curves are close one to another. At smaller radii, neutrino velocities are larger than CDM 
velocities. Note that Urms(?") of cold particles, which (to some extent) mimics the profile of the 
galaxies velocity dispersion in the clusters halo, stays almost constant, at least up to distance 
1.5h~^ Mpc. We define a /oca/ velocity dispersion crrms(?"), which mimics the local "temperature". 
Because there could be bulk fiows inside the clusters halo, some fraction of the peculiar velocities 
will not contribute to temperature. Let us define the bulk fiow velocity Vf[r) as the mean velocity 
of cold particles within a sphere of radius 150h~^ kpc centered on the each cold particle. Then 
the local velocity dispersion is CTrmsif) =< {v — u/(r))^ >^/^. The profile of CTrmsif) is shown as 
the dash-dotted curves in Figure 6. The difference between the two upper curves and lower curve 
Crms(?") is a measure of the departure from hydrostatic equilibrium in the cold and hot components. 
Figure 6 indicates that the hydrostatic equilibrium is a reasonable assumption within about radius 
~ lh~^ Mpc, but fails significantly at larger radii. It can be well explained by radial infall into the 
cluster halo. Also note that the cluster (a) is more distorted by the bulk fiows, associated with its 
close neighbor as compared to the cluster (b). 

5 Observational Tests on Clusters Density and Velocity 
Profiles 

In the CHDM model there are three matter components, cold, hot and baryonic, which are unequally 
distributed within the cluster. The baryonic component can be probed with the X-ray emission 
of intracluster gas. The distribution and motion of galaxies are primarily signatures of the cold 
component. The weak gravitational lensing would be produced by the total gravitational potential 
of the cluster. Therefore a combination of all three sorts of data would provide a test for the 
CHDM model. In this section we discuss some expected properties of clusters in the CHDM model 
regarding these observational tests. 



5.1 X-ray Surface Brightness and Temperature Profiles 

In this Section we use a simple model to check to which extend the dark matter density fit (13) 
is compatible with the X-ray data for clusters. By all means, this does not replace the numerical 
simulations with dark matter and gas, however, this is an instructive estimation. We will use, 
for simplicity, the model hydrostatic equilibrium of hot gas in the spherically symmetric static 
gravitational potential produced by the dark matter with the density profile (13). Then we will 
use well-known X-ray surface brightness profile for the cluster A2256, to calculate a theoretical 
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temperature profile. An actual temperature profile is compatible with a constant, the departure of 
the theoretical temperature from the constant measures the accuracy of this modelling. 



Let us assume that the gas density is spherically symmetric and is given by the /3-model (Cav- 
aliere & Fusco-Femiano 1976) 



Pgas 



1 + 



-3/3/2 



(14) 



where a is a core radius. Then the X-ray surface brightness, Ex(r), is 

-3/3-hl/2 

Ex(^) = Sx(0) ' 



1 + 



(15) 



where 6a is the angular size corresponding to a. The X-ray surface brightness of clusters typically is 
far from being spherically symmetric. However, the azimuthally averaged surface brightness can be 
fitted very well by formula (15). In Figure 7a we plot, as an example, Sx(^) for the Abell cluster 
A2256 (Briel & Henry 1994). This object is a rich, symmetric, partially relaxed cluster, similar to 
the Coma cluster. The best fit for the parameters fi and a are fi = 0.81 it 0.01, and a = 0.27h~^ 
Mpc. Note a very high accuracy ~ 1% in determining /3. 

Assuming gas is in steady state hydrostatic equilibrium with the spherically symmetric binding 
cluster gravitational potential, we have (Sarazin 1986) 



47r 



r Jo 



dr'r'^ p{r') 



kT[r) jdlnpgas{r) dlnT[r 



pnipG 



Jin r 



+ 



Jin r 



(16) 



where pnip is the mean particle mass in the gas, and T[r) is the gas temperature. This formula, 
in principle, relates the total dark matter mass density /o(r), given by (13), to the gas density 
(14). Unfortunately, even this oversimplified model based on (16) requires the knowledge of the 
temperature profile r(r), in order to reconstruct the dark matter density p{r). The measured X- 
ray temperature profile T[r) at present is known with a low accuracy, see Figure 7b. The actual 
temperature is compatible with the constant 7 Kev. The usual strategy is to insert the dark 
matter density in a chosen form (for instance, in the form (13)) into (16), and solve it for T[r) as 
function of the parameters of dark matter density profile. (Hughes 1989, Briel, Henry & Bohringer 
1992). One finds from (16) 



T{r) = pgUrr'{T{0)pgUO) 



ATTpnipG dr 



k 



r 



/2 



Pgas 



1 II 112 I ll\ 

dr r p[r ) 



:i7) 



The resulting T[r) profile is compared with the observed data to constrain parameters of the 
density profile. Despite the low accuracy of the measured temperature, this test turns out to 
be sensitive to the modelling of the dark matter profile. Our approach is complementary to the 
method of reconstruction of the X-ray surface brightness by adjusting the dark matter profile under 
the assumption that the actual temperature is constant, which is compatible with observations 
(Miralda-Fscude & Babul 1995). 
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Substituting the expressions (13) and (14) in (17), and using their asymptotes at large r it is 
easy to see that the integral in the right hand side converges to some constant J, which depends 
on i?, a,a,/3 and the dark matter density at the center po. The analysis shows that the integral is 
rapidly saturated at r ~ 0.5 Mpc, since the expression in the bracket goes rapidly to the constant 
/ = T[0)pgas{0) — J. Therefore at larger r the temperature T[r) crosses zero if / < 0, or increases 
very fast as p~}g{r) if / > 0. A physically acceptable temperature profile can be achieved only if 
r(0)/Jpas(0) is tuned to be close to J with an accuracy better than 1%. 

The projected temperature is 



where (Tgas{R) is the projection of the value p^gg_g{r). The theoretical temperature profile reduced 
from this formula for the chosen parameters of densities in gas (14) and total dark matter (13) 
components is plotted at Figure 7b. The data points correspond to the actual temperature mea- 
surement for the A2256 (Briel & Henry 1994). The theoretical prediction for T{r) from the CHDM 
model are compatible with the data for the A2256 cluster outside the radius ~ 0.5 Mpc, but depar- 
tures from the data up to factor 2 at the center. The temperature spike around the center depends 
on the details of the dark matter density peak. Traditional King profile for the dark matter density 
with the parameters adjusted to the numerical results of Fig. 4 gives much bigger departure of 
r(r) from the constant. Thus, either the model of hydrostatic equilibrium or the model (13) of the 
spherically symmetric dark matter distribution should be modified. 

5.2 Mass Density Surface Profile from Lensing 

The method, based on the weak gravitational lensing, can in principle, reconstruct the total pro- 
jected surface density E(r)tot. In the CHDM model, halos of galaxies are made from the cold dark 
matter particles. We expect therefore that in clusters, galaxies are tracers of the cold dark matter. 
From the known cold dark matter density profile, one can obtain the projected surface density of 
cold dark matter 



In the same way one can obtain the total projected surface density E(r)tot = Ti{r)cDM + S(r)j,. 

To see the difference in E(r)tot and Yj{r)cDM ^ In Figure 8 we plot Yj{r)cDM as dashed curve, 
and E(r)tot as the solid curve for the CHDM model with 2.3 eV neutrinos, calculated from formulae 
(19), (13) and /(r) from simple analytical model of Sec. 4.1. The lesson is that both total and CDM 
surface densities have almost similar shapes, the amplitude of the total surface density exceeds the 
CDM one by about 10%. The presence of the neutrino halo in the clusters basically changes only 
the amplitude of the surface density. Both uncertainties in the galaxies surface densities (Kaiser & 
Square 1993) as well as in the reconstructed surface density from the lensing (Fahlman et al 1994) 
make the neutrino halo of clusters in the fiat CHDM model with Oj, = 0.2 virtually undetectable 
by this method. 






(19) 
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However, in the framework of the CHDM model one can calculate the details of the mass density 
surface profile in the model. In Fig. 8 note the asymptotic slope E(i?)tot oc R~^'^ at large distance 
from the center, the asymptote E(i?)tot oc logR at small i?, and the bend of the surface profile at 
i? ~ 0.1 Mpc. Also note the remarkable difference between the hot gas mass density in Figure 7a 
and that of dark matter in Figure 8. In the ACDM model the density profile is steeper inside the 
cluster, therefore one expect the bend of the mass surface profile is at smaller R. The asymptotic 
slope —1.5 is the same. The clusters in the ACDM model are more relaxed than those in the CHDM 
model. Thus, the details of the surface mass density E(i?)tot might be another test to distinguish 
cosmological models. 



5.3 Cluster Galaxy Velocity Dispersion Profile 

The velocity dispersion cry[r) of galaxies as function of the radius has a significantly different profile 
in different models. As we described in Sec. 4, numerical simulations in the CHDM model show 
that cry[r) is almost constant in the tested range of scales up to distance ~ 3h~^ Mpc. The projected 
velocity dispersion is constructed from cry[r) by the formula similar to (18), and therefore also close 
to the constant. On the contrary, the ACDM model as well as the open CDM model show cry[r) 
decreasing with radius (c.f. Crone et al 1994). 

The observed projected galaxy velocity dispersion, apparently, has a central maximum, and 
then decreases with distance from the cluster center (e.g. Kent & Gunn 1982; Sharpes, FUis 
& Gray 1988). Potentially, this is very interesting test because this is in the confiict with the 
CHDM prediction, and but not with ACDM model. However, for more definite conclusion we 
have to take into account the complicated observational procedure of the rejecting of interlopes. 
The CNOC cluster redshift survey indicates a smaller fall of the projected velocity dispersion up 
to ~ 1.5h~^ Mpc (Yee et al 1994). (Note that this survey covers the range of clusters redshift 
0.2 < z < 0.5, where cosmological constant in the ACDM model has no significant cosmological 
impact). Thus, more data on the velocity dispersion profiles in clusters on different redshifts provide 
a discriminating test of cosmological model. 



6 Discussion 

We have studied several effects, responsible for cluster formation and structure in the Mixed Dark 
Matter models. Neutrinos in the range of masses up to a few eV contribute to the halo of clusters 
of galaxies, but not to that of galaxies or group of galaxies. Protoclusters form first from the cold 
dark matter particles, hot neutrinos are accreted afterwards onto the CDM protoclusters. Cold 
dark matter dominates most of the cluster mass density. The neutrino mass density is especially 
suppressed relative to the CDM density at the cluster center. At a distance of lh~^ Mpc the ratio 
of densities of both components saturates at the background ratio. Thus, only ~ 80% of dark 
matter is clustered inside the Abell radius. Therefore in the CHDM model the ratio of baryons to 
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the dark matter inside the clusters should be in 1.25 times greater than the background ratio. This 
is not enough to resolve the problem of large baryon fraction in clusters (Briel, Henry & Bohringer 
1992; White et al 1993). In the model with a cosmological constant, the ratio of baryons to the 
cold dark matter inside the clusters is about 3 times greater than the cosmological ratio of baryons 
to the total dark matter, including A-term. Thus, there is no baryonic catastrophe problem in the 
ACDM model. 

We found the density and velocity profiles of cold particles and neutrinos in clusters. The 
slopes a ~ 2.5 of the halo density profiles in both components are similar. The exact value of a 
depends on the neutrino mass. It should be noted that from the density profile slope alone one 
cannot distinguish the CHDM model, the ACDM model, or tilted CDM models. On the other 
hand, the velocity profile in the CHDM model is almost constant, while it is falling off with radius 
in the ACDM model. The projected velocity dispersion in clusters is a strong test to distinguish 
cosmological models. 

The oversimplified model of the spherically symmetric hot gas distribution in the hydrostatic 
equilibrium in the CHDM clusters was used to calculate the X-ray theoretical temperature T{r) 
from given X-ray surface brightness Sx(^) and spherically symmetric dark matter density profile. 
Around the center the theoretical temperature profile departures from the constant compatible 
with the data on A2256 cluster. 

The total mass density in the clusters outskirts has a ~ 20% excess relatively to that inside the 
cluster, due to the neutrino contribution. Galaxies in clusters, presumably, are good tracers of the 
cold dark matter component. The weak gravitational lensing would be a probe of the total cluster 
mass, due to cold dark matter and neutrinos. We found that the total projected mass density 
differs from the cold dark matter projected density by amplitude but not by the shape. The weak 
lensing better reconstructs the shape of the projected density. The slope of the total projected 
mass density of the clusters predicted by the model is ~ —1.5. The same slope is predicted in 
the ACDM model. The weak lensing method in principle, could probe the model by more subtle 
statistical tests of the cluster surface mass distribution. 

In summary, our study of the cluster properties in the CHDM model are adding several tests of 
the model. 

We thank S. Tremaine for useful discussions of the neutrino phase space density form. L.K. 
thanks the support from IFA, University of Hawaii, D.P was supported by CITA and CIAR cosmol- 
ogy program, and partly by the ITP UCSB. A.K. thanks IfA, University of Hawaii, for hospitality 
and the NSF grant AST-9319970. P.H. was supported by the NSF grant AST 9119216 and AST 
950515 and NASA grant NAG5-1880. 
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APPENDIX A: Linear Fluctuations in the CHDM model 



We have calculated the linear evolution of adiabatic perturbations in the model, treating neutri- 
nos kinetically through the Vlasov equation, and the CDM component as dust. We describe small 
scalar perturbations around a homogeneous, isotropic flat background universe in the form 



ds"^ = a{ri f [(1 + 2<!>)dri^ - (1 - 2^)dxadx"] 



(20) 



where a(?y) is a scalar factor of the universe, i] = / ^ is a conformal time, $ and ^ are gauge- 
invariant potentials describing scalar metric perturbations. In the homogeneous approximation, the 

-|- 1 ] , is the comoving momentum. 



neutrino distribution function is Fermi-Dirac Fo[q) - 
Let fl^{q",ri) be the linear perturbations of the neutrino Vi distribution functions. The Einstein 
equations for perturbations in Fourier decomposition (n„ is a comoving wave vector) are reduced 
to 



Stm a' \ I Sty a' 
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+ E 
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d q Jm'l.a'^ + + 
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.2^2 



d^'q 
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(21] 



Iml^a^ -\- q^ 

where ' 

The flrst-order corrections to distribution functions /^(<?",?7) satisfy the kinetic equations 



see Pogosyan & Starobinsky (1993) and Pogosyan & Starobinsky (1995b) for details. 



+ 



^ mla? -\- q^ 
IdFo 



J n 



q dq 



[22] 



Perturbations in the CDM (m) and radiative (r) components obey standard hydrodynamical equa- 
tions of motion. In terms of the invariant density perturbation 8e and velocity they are: 



'8t 



3 (8e^ 

4 \Z 



n 

+ — 
a 

a 



3*', (.r/a)'=2- + $. 



(23) 



These equations were solved numerically. 
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APPENDIX B: The Emden Equation for the CHDM Halo 



From the phase space density constraint, we have 

no < 0.5. (24) 
Substituting eq. (11) and (24) into (5), we obtain for the neutrinos mass density 



4 S 

m,,.a 



[l + exp (-//)) 



J 1 + exp (^^ - /ij 

Note, that actual density is smaller by factor 10. The parameter // can be fixed by considering the 
asymptote of the right hand side of eq. (25) at large radii. Indeed, at large r, from eq. (25) it 
follows that Pi,^{r) oc exp ( — ^(r)/(T^), similar to the CDM density (10). At large distances their 
ratio should matches the background ratio Oj,^ : VIcdm- From this we obtain 

l + e^^^,^lQ■-^■^-.pcDM{^)■[-^] . (26) 

Thus, in our model the integral (25) mildly depends on the CDM central density. Finally, the 
neutrino density profile can be estimated by formula 

Oj,. /2 / X /■°° dyy'^ 

ftcDM^ Jo exp/i + exp(j/V2 + (^(r)/(72) 

where p is given by equation (26). 

Poisson's equation for the gravitational potential (f>[r) is 

+ = iPcDM{r) + p,{r)) , (28) 

where the CDM and neutrino densities as functions of cf){r) are given by eqs. (10), (27). It is 
convenient to introduce the dimensionless variables, for the potential 

<f>{r) = -u{r)(j\ (29) 

and for the length 



zcr 



^iirGpcDMiO) 

Poisson's equation is reduced to the form 



(30) 



^ 2du ^ ^ 0., 2 ro- dyjl 

dz^ z dz flcDM\T^Jo exp p -\- exp (y^ /2 — u{z)) 

The boundary conditions at the center z = are u = and du/dz = This is a modification of 
the standard form of the Fmden equation. In the limit of large r the last equation has the solution 

PcDM{r)^-^Jl + ^^). (32) 
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Our toy analytic model is based on simple phase space density distributions (9) and (11), which 
lead to the mass density distributions (10) and (27), similar to the isothermal profiles, see Fig. 
3. This allows us to understand an important feature: the central mass density of neutrinos is 
significantly suppressed by the factor / relatively to the CDM density. Our numerical "coarse- 
grained" mass density has similar properties, see sec. 4. However, as numerical simulation of sec. 
4.2 shows, the realistic cold component density profile (13) is quite different from that of Fig. 3. 
In principle, the actual central CDM mass density might be very high. There is an interesting 
question if the central mass density of the hot component is suppressed by factor / but also high? 
Unfortunately, our analytic model as well as numerical simulation do not address this question. The 
finding of the phase space density distribution which could correspond to the mass density (13) is 
much more complicated problem, and numerical study requires much higher resolution. However, 
we believe that the neutrino central mass density has no cusp and cannot be very high. This is 
because the neutrinos begin clustering at relatively late stages, well after the clustering of the cold 
component. 
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Figure Captions 



Figure 1: The transfer function C{k) of the total density fluctuations in the CHDM model 
with Oj, = 0.2. An upper curve corresponds to the models with only one 4.6 eV massive neutrino; 
a middle curve - two 2.3 eV neutrinos; and a lower curve - three 1.53 eV neutrinos. 

Figure 2: Typical density profiles around the peak, solid curve for CDM and dashed curve for 
2.3 eV neutrinos. The upper panels corresponds to the top-hat filtering of the initial density field 
on the scale 1O^"^M0; the peak height is p = 1.68. The lower panels corresponds to the filtering on 
the scale lO^^M©, the peak height p = 2.5. The left panels corresponds to the redshift z = A; the 
right - z = 0. 

Figure 3: The left panel shows the ratio of neutrino and CDM density profiles /(r) = 
Pi,{r) I pcdm{''^) as function of r. The solid curve corresponds to 2.3 eV neutrinos; dotted curve 
to 4.6 eV neutrinos, and dashed curve to 1.56 eV neutrinos. The right panel shows three pairs of 
cold (upper curves) and hot (lower curves) particles density profiles for these three models. 

Figure 4: The left panel shows the density profiles for the cluster (a) with the mass 6.5/i~^ lO^'^M©; 
the right panel shows the density profiles for the cluster (b) with mass 5.5/i~^lO^'*M0. The upper 
panels show the ratio of the cold and hot particles densities; normalized in such a way that unity 
corresponds to the background ratio. The lower panels shows CDM density profiles and parameters 
of the fitting form (13). The distances are given for h = 0.5. 

Figure 5: The same as in Figure 4 but for the "old" CHDM model with one 7 eV neutrinos. 

Figure 6: The profile of the velocity dispersions Urms(?") in the clusters. Left panel is for cluster 
(a), right panel for cluster (b); both for the 2.3 eV model, see Figure 4. The solid curve shows 
Vrmsi''^) of cold particlcs; the dashed curves that of the hot particles. The local velocity dispersion, 
(Trms{i^) =< {v — u/(r))^ >"^''^, is showu as the dash-dotted curves. 

Figure 7: (a) Observed X-ray surface brightness Sx(^) for the Abell cluster A2256 (Briel & 
Henry 1994); solid curve shows the /3— model fit (15). (b)The measured X-ray projected temperature 
r(r) for this cluster. Solid curve is the projected temperature obtained from (17). 

Figure 8: Solid curve is the total projected surface density E(r)tot (in units of mean density 
by Mpc); dashed curve is projected surface density of the cold component Yj{r)cDM only. 
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